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Abstract 

The paper Is a review of sound generation in 
turbulent shear flows. It emphasizes simultaneous 
calculation of the turbulent flow along with the 
resulting sound generation rather than the alterna- 
tive acoustic analogy approach. The first part of 
the paper Is concerned with solid surface Interac- 
tion. The second part concentrates on the sound 
generated by turbulence Interacting with Itself. 

Introduction 

The subject of this paper became a serious 
scientific discipline about 35 years ago when 
Llghthill 26 * 27 published his acoustic analogy 
theory of jet noise. That work has more or less 
dominated the subsequent development of this sub- 
ject, which Is still somewhat Incomplete— even 
though It has undergone little change In the past 
five years or so. This paper does not place much 
emphasis on the acoustic analogy but rather concen- 
trates on an alternative approach, which may be 
more readily adapted for use on large-scale com- 
puters to obtain more detailed Information about 
the sound field than would be possible from the 
acoustic analogy. 

This approach amounts to little more than cal- 
culating the unsteady flow that produces the sound 
simultaneously with the resulting sound field, 
starting from some prescribed upstream state that 
is Ideally specified Just ahead of this region where 
the sound generation takes place. To make progress 
without resorting to full-scale numerical computa- 
tion requires that the governing equations be line- 
arized about some appropriate mean flow. But that 
ultimately has to be done, either Implicitly or 
explicitly, even with the acoustic analogy approach. 
I have no doubt that the day will come when 
turbulence-generated sound Is calculated directly 
from the Navler-Stokes equation, but, to my know- 
ledge, that has yet to be done, and there Is much 
to be learned from the existing work, which has 
often led to relatively simple formulas that show 
encouraging agreement with experiment and produce a 
great deal of Insight Into and physical understand- 
ing of the sound generation process. This paper Is 
a very selective review of that work, and It adopts 
a very specific viewpoint. 

The use of linearized theory to calculate tur- 
bulent flows or, better yet, changes In turbulent 
flows Is a branch of turbulence theory now known as 
"rapid-distortion theory" (see Moffatt 31 ). It 
assumes that the following assumptions are satisfied 
(Hunt 24 ): (a) that u'/U « 1, where u' Is the rms 
turbulence velocity and U Is the local meanflow 
velocity; and (b) that the Interaction or change 
being calculated be completed In a time, tj say, 
that Is short compared with x<j e cay» where T<j e cay 
Is the decay time or lifetime of a typical turbulent 
eddy 0(l/u')~ l being the characteristic size of 
turbulent eddies. Rapid-distortion calculations are 
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usually based on the Invlscld equations-- an approx- 
imation that Is justified when both the mean flow 
and turbulence Reynolds numbers are large. The 
Important point here Is that the radiated sound 
field can be determined as a by-product of any such 
rapld-dlstortlon-theorv calculation, as long as 
compressibility effects are retained . 

While It might seem most logical to begin by 
omitting solid surface effects and to Include them 
only after the turbulence self-noise problem has 
been appropriately dealt with, it turns out that 
the solid boundaries actually simplify the problem 
and allow a more rigorous treatment In, at least, 
some cases. Consider then a hlgh-Reynolds-number 
turbulent alrjet such as shown schematically In 
Fig. 1. The maximum turbulence level occurs along 
the centerline of the Initial mixing layer, Indica- 
ted by the dashed line In the figure. Here the 
ratio of the rms turbulence velocity to the local 
mean-flow velocity Is roughly equal to 0.24 
(Bradshaw et al. 5 ) which, although not all that 
small, would probably still be considered to be an 
acceptable "small parameter" to many classical 
applied mathematicians. Condition (a) Is therefore 
reasonably well satisfied. 

Solid Surface Effects 

Now suppose that a semi-infinite but Infinites- 
imally thin flat plate Is Inserted Into the flow as 
shown In the figure. Then the Interaction between 
the turbulence and the leading edge will be comple- 
ted In a time tj = 0(l/U), which. In view of the 
smallness of the turbulence Intensity, Is fairly 
small compared with T^ecay * 0(l/u'). Thus, Invls- 
cld rapid-distortion theory applies, and the Inter- 
action between the turbulence and the edge can be 
calculated by linearizing the Invlscld equations 
(the Euler equations) about the mean flow. 

Since the ratio of the cross-stream to stream- 
wise components of the mean-flow velocity Is of the 
order of (u'/U) 2 (Tennekes and Lumley 40 ), the order 
of approximation will certainly not be diminished 
If this flow Is taken to be a unidirectional trans- 
versely sheared flow. The Important advantage of 
using this flow Is that It Is Itself a solution of 
the Invlscld equations (for any velocity profile). 
The resulting expansion Is then a rational pertur- 
bation that can, In principle, be carried to arbi- 
trary order without Internal Inconsistency. The 
lowest order equations are the same as those used 
In Invlscld stability theory, l.e. the Rayleigh 
equations (see Betchov and Crlmlnale 4 ), and as 
already Indicated, the radiated sound field can be 
determined as part of the solution to these 
equatlons^-provlded, of course, that compressibility 
effects are retained. 

To accomplish this, one must first decide on 
an appropriate representation of the Incident turbu- 
lence. This would be rather easy to do If the mean 
flow were completely uniform, since any solution 
for the unsteady velocity/pressure fluctuations 
could then be decomposed Into the sum of an 


“acoustic solution" that carries no vortlclty and a 
vortical solution which produces no pressure fluc- 
tuations and which Is often referred to as the 
"gust" or “hydrodynamic" solution. The latter Is 
used to represent the Incident turbulence In most 
problems that Involve the Interaction of turbulence 
with solid surfaces embedded In uniform mean flow. 
Its suitability for this purpose Is largely due to 
the following reasons: 

(1) It does not become Infinite anywhere In 
space, even In the absence of solid surfaces, so 
that it can describe the turbulence field that would 
exist If the surfaces were not present. 

(2) It Involves two arbitrary "convected" quan- 
tities that can be specified as upstream boundary 
conditions to describe the turbulence entering the 
interaction zone In any given problem. This seems 
to be the appropriate degree of generality, because 
the the vortlclty Is a convected quantity that has 
only two Independent components (since Its diver- 
gence must vanish) . 

(3) It has no acoustic radiation field at sub- 
sonic speeds and will In fact vanish exponentially 
fast at transverse Infinity If the mean and unsteady 
vortlclty fields are sufficiently compact. 

Decomposition of the solution Into acoustic 
and vortical parts Is no longer possible when the 
mean flow Is nonuniform, but the compressible 
Rayleigh equations still possess a solution that 
has the three properties listed above and. In fact, 
approaches the "vortical" solution on a uniform 
mean flow In the limit as the mean flow approaches 
a uniform flow (Goldstein 18 ' 19 and Mohrlng 32 ). This 
would then seem to be the natural generalization of 
the latter to nonuniform flows, and It would there- 
fore seem appropriate to refer to It as the "gust" 
or "hydrodynamic" solution and, more Importantly, 
to use It to represent the Incidence turbulence. 

We therefore use It to represent the Incident 
turbulence In the present problem of a large flat 
plate embedded In a turbulent shear flow. Since 
the gust solution does not satisfy the zero normal 
velocity boundary condition at the plate. It Is 
necessary to add another solution to cancel this 
component of velocity. Unlike the gust, this latter 
solution does not vanish exponentially fast at 
Infinity but rather behaves like a propagating 
acoustic wave there (Goldstein 19 ). In other words, 
the plate Is able to "scatter" the nonpropagating 
motion associated with the gust Into a propagating 
acoustic wave (Ffowcs Williams and Hall 13 ). 

The problem also possesses an elgensolutlon 
associated with the spatially growing Instability 
wave that can propagate downstream from the edge on 
the Inflectional mean-velocity profile (Crlghton 
and Lepplngton 9 and Goldstein 20 ). The solution Is 
therefore not unique! It could be made unique if we 
required that It remain bounded at Infinity (since 
that would eliminate the elgensolutlon which grows 
without bound there). But since the linearization 
Is only valid In the vicinity of the leading edge. 

It Is probably not appropriate to Impose a 
"boundary" condition far downstream In the flow 
where all sorts of nonlinear effects will have had a 
chance to Intervene (Rlenstra 36 ). 


One can therefore look for an alternative way 
to make the solution unique. This can be done by 
treating the steady-state solution, which Is, of 
course, the one of Interest here, as the long time 
limit of the solution to an Initial-value problem, 
and then Imposing a "causality condition" In the 
sense that the solution Is required to be Identi- 
cally zero before the Initial time when the Incident 
disturbance Is "turned on" (Crlghton and 
Lepplngton 9 ) . 

But Rlenstra 36 argued that an Initial condition 
Imposed In the distant past may not be relevant to 
the steady-state solution, since the linearization 
might only be valid over a relatively short Interval 
of time. One might therefore consider a third way 
of making the solution unique. This amounts to 
using the elgensolutlon to eliminate the leading- 
edge singularity that appears In both the bounded 
and causal solutions, l.e. by satisfying a leading- 
edge "Kutta condition" (Goldstein 20 ). This may 
be rationalized by noting that the Instability wave 
represents downstream vortex shedding that could 
adjust Itself to eliminate the singularity In the 
Invlscld solution and thereby prevent any flow sep- 
aration that would otherwise occur at a very sharp 
edge. 

It Is not entirely clear which of these three 
solutions Is correct, but I expect Rienstra's 36 
argument Is Invalid and that Imposition of causality 
Is probably appropriate here. Goldstein 19 compared 
the theory with the data of Olsen, 35 who measured 
the sound radiated In one-third-octave frequency 
bands as a function of the angle from the jet axis 
In a plane perpendicular to that of the plate. 
Comparison of experiment and theory Is shown in 
Fig. 2. The top of the figure corresponds to the 
high-frequency limit where the Instability waves 
are "cut off" and the Issues of causality and Kutta 
condition are Irrelevant. However, the low- 
frequency causal solution, which Is shown at the 
bottom. Is strongly affected by the Instability 
wave. The agreement between experiment and theory 
Is good, but the causal and leading-edge Kutta con- 
dition solutions have the same low-frequency limit, 
and one cannot conclude from this comparison which 
Is correct. However, the bounded solution behaves 
quite differently In this limit and consequently 
does not agree with the data. 

Sound Generated by Turbulence Interacting 
With Itself: The Jet-Noise Problem 

Having achieved some success In using linear 
theory for the turbulence-leading-edge Interaction, 
It Is natural to try using It to calculate the sound 
generated by turbulence Interacting with Itself, 
l.e. to deal with the problem of jet noise. I have 
already pointed out that the ratio rms turbulence 
velocity to local mean-flow velocity Is reasonably 
small In the region of maximum turbulence level, so 
that the first requirement for the validity of 
rapid-distortion theory (see Introduction) Is satis- 
fied. However, the Interaction time tj, which In 
the present context should be taken as the time for 
the sound generation to occur. Is now equal to the 
decay time T<j e cay of the turbulence, and the 
second requirement Is not. But with no better 
alternative at hand, we might still attempt to 
Introduce the same small parameter as before, l.e., 
u'/U, and carry the corresponding asymptotic expan- 
sion to Its logical conclusion. Like the more ad 


2 



hoc acoustic-analogy approach, this systematic pro- 
cedure assures that all appropriate conservation 
laws will be satisfied and that the acoustic sources 
will be of the appropriate multipole order. But It 
seems to have certain advantages over the acoustic 
analogy In that It provides a "rational" framework 
for assessing the Internal consistency of the vari- 
ous jet-noise analyses. It may also apply to some 
physically realizable flow, which Is hopefully not 
too different from the real turbulent flow of Inter- 
est, and finally. It provides a method for Identify- 
ing acoustic sources and distinguishing acoustic and 
nonacoustic components of the unsteady motion. 

The Basic Equation 

The lowest-order equations are, on the face of 
It, the same as before, l.e. they are the compress- 
ible Rayleigh's equations. It Is well known 
(Betchov and Crlmlnale 4 ) that the velocity compo- 
nents can be eliminated between these equations to 
obtain a single equation for the normalized first- 
order pressure fluctuation = Pi/p 0 C§. where 
Is the actual first-order pressure fluctuation, 
PoflTt) the mean-flow density, and C 0 (lT t ) Is the 
mean-flow sound speed where the latter two quanti- 
ties depend only the cross-stream coordinate jcj * 
(X 2 . X 3 } - X], X 2 , X 3 denoting Cartesian coordi- 
nates with In the mean flow direction. This 

equation can be written symbolically as 

Ln*i * o f (l) 

where L denotes the third-order linear-wave 
operator 

L = St(^2- V - C S*) +2C 0< 7 U) "1^. (2) 

t denotes the time, and D/Dt = 3/3t+U3/3xi Is the 
convective derivative based on the mean-flow 
velocity UCxt). 

Since solid boundaries are acoustically irrele- 
vant for the turbulence self-noise problem. It Is 
appropriate to suppose that the flow Is defined 
over all space. Then {for reasons given In the 
section on solid surface Interactions) the "gust” 
or "hydrodynamic" solution Is an appropriate solu- 
tion of Eq. (1). But this equation also has (spa- 
tially growing) Instability-wave solutions which 
can exist whenever the mean flow Is Inflectional 
(Betchov and Crlmlnale 4 ). Since many Investigators 
have argued (Crlghton , 8 * 7 Liu , 29 Tam and Chen , 38 
Haertlg , 22 Gaster, Kit, Wygnanskl , 44 and others) 
that these latter solutions correspond to the exper- 
imentally observed large-scale turbulent structures. 
It would seem appropriate to Identify the gust solu- 
tion with the "fine-grained" (or relatively "fine- 
grained") turbulent motions. 


carried to the next order If It Is to be used to 
calculate radiated sound. The second-order normal- 
ized pressure fluctuation again satisfies a 
third-order wave equation, but It Is more convenient 
to work with the Isentroplc density fluctuation 

k - 1 2 

n = n 2 - — n r ( 3 ) 

where k Is the specific heat ratio. Then n 
satisfies 

Ln = Y* (4) 

which, except for the Inhomogeneous source term 

* = ot v * ? - 2 fr/ vu - < 5 > 

Is the same as Eq. (1) for the first-order normal- 
ized pressure fluctuation. 

Equation (5) Is Identical to the source term 
that would be produced by an externally applied 
fluctuating force per unit mass T = {^^ 2 ^ 3 } and 
might therefore be thought of as a dipole-type 
source, since a fluctuating force produces such_a 
source when there Is no mean flow. The force T 
Is not, of course, arbitrary but Is now given as a 
quadratic function of the first-order solutions, 
viz. 

< 1.3 - 

(6) 

where u(^ are the first-order velocity 
1 2 

fluctuations and = k^?T^ Is the first-order 

square sound speed fluctuation (# being the gas 
constant and T-j being the first-order temperature 
fluctuation) . 

Equations (4) to ( 6 ), with some relatively 
minor differences, were first derived by Lllley , 28 
who used the acoustic-analogy approach. The result 
Is now commonly referred to as Lllley's equation. 

In the present approach, It arises as the equation 
for the composite second-order pressure fluctuation 
n with a source term y that Involves only first- 
order solutions. Since these solutions, which sat- 
isfy the homogeneous Eq. (1), have no acoustic 
fields at subsonic speeds, while the second-order 
solution does, our expansion provides a (conceptual 
If not experimental) mechanism for Identifying 
acoustic and nonacoustic parts of the unsteady 
mot 1 on--but there are some complications. 

The Sources of Sound 


However, there are experimentally observed, 
large-scale motions which, on a global basis, seem 
to bear little resemblance to any motion that can 
be represented by either the gust or linear 
Instability-wave solutions. This should come as no 
surprise, since we have already noted that the lin- 
earized solution can at best remain valid over rela- 
tively small streamwlse distances. 

We have seen that the gust solution produces 
no acoustic radiation at subsonic speeds, and the 
same can be said for the Instability waves (but see 
below). The asymptotic expansion must therefore be 


The second term In Eq. ( 6 ) represents a dipole- 
typen source due to the temperature fluctuations In 
the flow (Tester and Morfey 42 ). While this source 
Is of real significance In actual high-temperature 
jet exhausts, I will not discuss it in this paper. 

I concentrate rather on the first term, which, being 
the divergence of the fluctuating (first-order) 
Reynolds 


stress uj 1 ^. corresponds to the source that 

would be produced by an externally applied fluctua- 
ting stress field. It might therefore be 
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Interpreted, by analogy with the zero-mean-flow 
case, as a quadrupole-type source. 

This latter term can be further decomposed Into 
a number of subsources by separating the first-order 

solution Into Its gust and linear (spatially 

growing) instability-wave components and, as before, 
Identifying the gust with the fine-grained turbulent 
motion. Unfortunately, this procedure cannot be 
carried to Its logical conclusion because the linear 
instability waves, which grow without bound on a 
parallel mean flow, will ultimately produce an 
unbounded source term In Eq. (4). It would then be 
Inappropriate to use this equation to calculate the 
acoustic field, since It Is Its global, and not Its 
local, solutions that must be used In such a calcu- 
lation. However, the real flow Is only locally 
parallel, and the slowly varying (rather than the 
parallel-flow) approximation should be used to 
represent the Instability waves, as was done by 
Crlghton and Gaster, 8 Tam and Morris, 39 and others. 
Then the source term In Eq. (4) will remain bounded, 
since the local growth rate of the Instability wave 
varies with the thickness of the jet or shear layer 
(It first Increases, reaches a maximum, and then 
becomes negative as the thickness Increases). 

However, supersonically traveling waves can be 
produced as a by-product of this approximation, and 
these latter waves will couple to the radiation 
field (Tam and Morris 39 ) when the first-order solu- 
tion Is rendered uniformly valid (by using an appro- 
priate singular perturbation procedure, such as the 
method of multiple scales 34 ). Our previous comment 
that the first-order solution has no radiation field 
therefore needs to be qualified. We return to this 
below, but for now the Important point Is that the 
first term In T should then describe the sound 
generation due to the following types of 
Interactions: 

(1) Linear Instabl 11 ty wave-fine-grained 
turbulence. 

(2) Linear Instability wave-linear Instability 
wave. 

(3) Fine-grained turbulence-fine-grained 
turbulence. 

This list may be Incomplete, of course, or 
even Inappropriate, since as I already Indicated, 
there are other types of large-scale motions In the 
jet that do not seem to be globally representable 
by either the gust or Instability wave. 

In any case. It Is clear that this list should 
only be taken as an Indication of the types of 
Interactions that can occur and should not be 
considered to be the result of a rigorous analysis. 
In fact we shall eventually show that nonuniform- 
ities In the asymptotic expansion cause these Inter- 
actions to occur at different asymptotic orders than 
the present formal expansion would suggest. One 
might then choose to Ignore the list entirely and 
argue that the experimentally observed turbulent 
motions should be used In place of the first-order 
solutions that appear In the source term (5), which 
Is, In effect, what Is done In the acoustic-analogy 
approach. However, I do not think that It should 
be dismissed entirely and therefore consider It In 
some detal 1 . 


Its first Item, l.e. the Instability wave fine- 
grained turbulence Interaction, has only been con- 
sidered very briefly In the literature (Ffowcs 
Williams and Kempton, 14 Berman, 3 Liu 29 ) and, to myl 
knowledge, only limited quantitative results have 
been obtained for this Interaction. It can be 
thought of as the sound generated by the fine- 
grained turbulence shaking the Instability waves 
and Is likely to emerge as an Important source mech- 
anism In relatively low Reynolds-number flows. 

The Instability wave-lnstabl 1 Ity-wave Inter- 
action may be related to the vortex-pairing events 
that occur In the Initial mixing region of a high- 
speed jet. These events can be experimentally 
enhanced by exciting the jet with an external acous- 
tic source tuned 1 to the most unstable frequency of 
the shear layer at the nozzle lip, as was done by 
Klbens. 25 He found this to cause a supression of 
the natural broadband noise of the jet with most of 
the sound being generated at subharmonics of the 
excitation frequency. By taking measurements In 
the near and far fields, Klbens 25 showed that there 
was no Doppler shift In frequency. Indicating that 
the sound was generated by nonconvectlng sources 
within the jet, whose locations he subsequently 
Identified with the vortex-pairing locations. 

However, quadratic Interactions between two- 
dimensional (or between axlsymmetrlc) Instability 
waves produce only subsonlcally traveling waves on 
a subsonic parallel mean flow, and these waves do 
not radiate sound. But the straightforward pertur- 
bation analysis of these Interactions is (like the 
straightforward nonparallel mean-flow analysis) 
nonunlformly valid In the streamwlse direction 
(leading to the so-called Kelley 45 resonance), 
and (as In the nonparallel analysis) supersonically 
traveling waves are produced when the straight- 
forward asymptotic solution Is rendered uniformly 
valid In that direction. The sound field can then 
be calculated by using a procedure similar to the 
one used by Huerre and Crlghton 23 for the sound 
generated by the nonlinear saturation of a single 
Instability wave. But a more systematic approach 
might be to adapt the procedure used by Tam and 
Morris 39 to this case. 

The Llghthl 1 1 Result 

The fine-grained turbulence-fine-grained- 
turbulence Interaction Is essentially the mechanism 
originally considered by Llghthl 1 1 . 2 °* 27 Difficul- 
ties such as those discussed In conjunction with 
Instability wave-lnstablll ty-wave Interaction may 
also occur when the present perturbation approach 
Is applied to this case. Llghthl 1 1 • s acoustic 
analogy theory leads to a stationary medium (l.e. 
classical) wave equation. He suggests that It 
should be possible to neglect variations In retarded 
time across the turbulent eddies (or correlation 
volumes) In this case, since the time l/C 0 
(1 - M c cos e) for a sound wave to cross a turbulent 
eddy will be small compared with the characteristic 
time l/u' of the sound source at subsonic jet 
velocities (Goldstein 17 ). Here M c Is the “convec- 
tion Mach number*' of the turbulence, and e Is the 
angle between the downstream jet axis and the line 
connecting the source point and the observation 
point. 

It Is a consequence of the linearity of the 
wave operator on the left side of Eq. (4) that the 
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sound radiated by any given turbulent eddy will 
then be Independent of that radiated by any other 
eddy. Llghthlll therefore argued that each eddy 
should behave like a point quadrupole source moving 
downstream with the "convection velocity" of the 
turbulence and that the entire sound field of the 
Jet could then be estimated by calculating the sound 
radiated by a "typical" turbulent eddy. This pic- 
ture turns out to be a slight oversimplification 
and was later corrected by Ffowcs Williams. 12 

In the L1ghth11! 26 /FfQwcs Williams 12 result, 
the mean-square pressure p 2 radiated In any pro- 
portional frequency band at a fixed source frequency 
Q, where u * Q(1 - M c cos ©) Is the actual fre- 
quency of the sound, behaves like 



(7) 


so that Its "directivity pattern" Is primarily 
determined by the Doppler factor (1 - M c cos 9) 
raised to the -5 power. These five Inverse Doppler 
factors produce a highly directional radiation pat- 
tern at high subsonic Mach numbers--wh1ch Is remark- 
ably similar to experimental observation. 


Solutions of Llllev's Equation 


Solutions of Eq. (4) with y treated as a 
moving point source can be Interpreted as correc- 
tions to the L1ghth111 26 /Ffowcs Williams 12 result 
(7) that account for the effects of the nonuniform 
surrounding mean flow. A number of us (Manl, 30 
Balsa, 1 * 2 Berman. 3 Goldstein, 15 * 16 * 21 Tester and 
Morfey, 42 Scott, 37 Lllley, 28 etc.) therefore decided 
to calculate the acoustic radiation from point quad- 
rupole sources moving through transversely sheared 
mean flows. The relevant solutions usually had to 
be obtained numerically, but relatively simple 
closed-form (or nearly closed-form) solutions were 
obtained In the lowand high-frequency limits uD/Uj 
« 1 and wD/Uj » 1, respectively, where D 
denotes the jet diameter (see Fig. 1) and Uj 
denotes the jet velocity. 


Low-frequency solutions were obtained for a 
round jet with arbitrary mean-velocity profile by 
Goldstein 15 * 16 and Balsa. 2 All components of an 
Idealized quadrupole convectlng through a stationary 
medium exhibit directivity patterns given by Inverse 
Doppler factors times sines and cosines of the 
observation angle. The low-frequency analyses show 
that only the x-|-x-j and x-|-r quadrupole compo- 
nents (where r Is the radial coordinate) retain 
this property In the presence of a parallel but 
nonuniform mean flow. The remaining quadrupole 
components exhibit directivity patterns given by 
more-complex formulas Involving the complete mean- 
velocity profile and the location of the sources 
within the jet (Goldstein 16 ). 

However, the low-frequency analyses uncovered 
the very surprising result that the mean flow causes 
certain quadrupole components to emit sound much 
more efficiently than they otherwise would— the 
mean-square pressures In the absence and presence of 
the mean flow being respectively, 0(Q 4 ) and 0(Q 2 ) 
as Q -> 0. The acoustic field of the x-|-r quad- 
rupole, which Is the only one of these more effi- 
cient sources that can be expressed In simple 
Ooppler-factor form. Is proportional to the local 


mean-velocity gradient. It Is worth noting that 
this source arises as much from the first member of 
the source term (5) as from the second, even though 
the former does not explicitly Involve the mean- 
velocity gradient (Balsa 2 ). 

Observed low-frequency jet-noise directivity 
patterns therefore depend on complex properties of 
the jet turbulence and mean flow that are difficult 
to estimate. But the mathematical results are con- 
sistent with the experimentally observed result 
that they will always be more directional than 
Llghthlll's Inverse five Doppler factors would 
Indicate, or, to be more specific, the analytical 
and experimental results show that the low-frequency 
sound should be more concentrated on the down-stream 
axis than Llghthlll's result Implies, with the on- 
axis sound being produced by the quadrupoles with 
one axis In the streamwlse direction. 

The high-frequency solutions, which were i 
obtained by Lllley, 28 Tester and Burrin, 41 balsa,'* 4 - 
Berman, 3 and Goldstein, 11 * 21 exhibit a "zone of 
silence" on the downstream jet axis. The acoustic 
field Is exponentially small In that region, which 
Is circumferentially asymmetric when the jet Is 
nonaxlsymmetrlc and/or the sound source Is located 
off-axis. It will fill the entire range of circum- 
ferential angles when e Is sufficiently close to 
the downstream jet axis (see Fig. 3), but will only 
occupy a limited range of angles (say <f>Cm1n < v < 
<Pmax> at larger values of e (say e Cm1n 
c © < ©cnax)* and finally It will disappear com- 
pletely when © > ©cmax* (See 4 «) These 
remarks only apply to subsonic Isothermal jets with 
monotonic or nearly monotonic mean-velocity pro- 
files. A host of complex Interference effects may 
occur when these restrictions are relaxed. 


As In geometric optics, the sound propagates 
along distinct rays In the high-frequency approxi- 
mation. Only one ray can reach the observer when 


© > ©cmax* but there w111 be at "least two rays 
reaching the observer when ©cmln < e < ©Cmax — a 
direct ray and a ray reflected from the boundary of 
the zone of silence (Goldstein 21 ). The correspond- 
ing sound waves can then Interfere, but the Interfer- 
ence term will be a rapidly oscillating function of 
angle and, since all acoustic measurements Involve 
some for of a spatial averaging, may not be experi- 
mentally observable. 


The mean square pressure radiated In any pro- 
portional frequency band of fixed source frequency 
Q will then be the sum of the mean-square pressures 
for each ray reaching the observer. The result for 
a convectlng quadrupole source, corrected for the 
Ffowcs Williams 12 effect. Is given by (Goldstein 21 ) 

3 2 

_ p - ft5 

p 2 a Lilji 1 - A(q>), 

(4irRC Qo C o ) <C (l - M COS ©r(l - M c COS ©r 

( 8 ) 

where R Is the distance between the source point 
and the observation point shown In Fig. 4, Pa3 and 
C« are the density and sound speeds at Infinity, 
and M Is the Mach number based on the mean flow at 
the source location and the speed of sound at 
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Infinity. The Q^j 
rupole strengths, 

denote 

the 

relative quad- 


v-j = q Q cos X, 

= q o sm 

X, 

\> 3 = cos e, 

(9) 

J (— 

> 

M cos e 

\ 2 

2 

cos e, 

(10) 


c /c > 

0 / 

1 - 


X Is the Initial circumferential angle made by the 
acoustic ray associated with Eq. (8), and A 
denotes a "circumferential directivity factor" that 
depends on the circumferential observation angle <p, 
the location of the sound source within the jet, and 
the mean-velocity and temperatures profiles of the 
jet. 

Equation (8) Is an exact high-frequency result 
that applies to jets of any cross section and with 
any transverse mean-velocity and temperature pro- 
files, but x and A must be calculated by solving 
a second-order ordinary differential equation In 
the general case. They are, however, given by rel- 
atively simple analytic formulas (Goldstein 2 ') for 
off-axis sources at arbitrary locations In a cir- 
cular jet with arbitrary velocity and temperature 
profiles . 

The circumferential directivity factor can be 
used to study the effect of nonaxlsymmetrlc jet 
velocity and temperature profiles In reducing jet 
noise below the flight path of a jet aircraft, which 
Is of considerable Interest for technological appli- 
cation (Von Glahn and Goodykoontz 43 ) . But for 
the present purpose, It Is appropriate to concen- 
trate on the azimuthal directivity pattern, which 
Is relatively unaffected by this factor. 

Equation (8) shows that the Inverse-Ooppler- 
factor exponent Is Increased from 5 to 7 In the 
high-frequency limit, since the local mean-flow 
Mach number and the turbulence-convection Mach num- 
ber are usually not very different In the regions 
of peak turbulence Intensity. This taken by Itself 
would cause the high-frequency sound (like the low- 
frequency sound) to be more directional than 
Llghthlll's Eq. (7) would predict. If, however, 
following Balsa, 2 the quadrupole Is assumed to be 
Isotropic so that 


of course, highly unlikely that the actual quadru- 
poles will precisely satisfy Eq. (11), but the anal- 
ysis strongly suggests that quadrupole component 
Interference effects can greatly reduce jet-noise 
directivity. 

Sound Generation Due to Streamwlse 
Variations In Mean Flow 

The formal asymptotic expansion In powers of 
(u'/U) can be continued to the third order. At 
this stage. Interactions between the first-order 
perturbation solution and the streamwlse variations 
In the mean flow will appear In the source term. 
[Recall that the ratio of the cross-stream to 
streamwlse components of the mean-flow velocity Is 
0((u'/ll) 2 ), while the first-order solution Is 
O(uVU).] Then, by decomposing the first-order 
solution Into Its gust and Instability-wave compo- 
nents and making the connection between the gust 
and fine-grained turbulence that we discussed above, 
we Infer that the source term now describes sound 
generation due to (a) the fine-grained turbulence 
Interacting with streamwlse variations (l.e. spread- 
ing) of the mean flow and (b) the Instability waves 
“Interacting with the streamwlse mean-flow varia- 
tions. The first mechanism has, to my knowledge, 
not yet been considered In the literature. The 
second will already be accounted for In the first- 
order analysis If the slowly varying approximation 
Is used to describe the Instability waves. This 
mechanism has been analyzed In a more ad hoc fashion 
by Crow and Champagne, ^ Ffowcs Williams and 
Kempton, 14 Huerre and Crlghton, 23 and Liu, 29 and In 
a systematic way by Tam and Morris. 39 However, Tam 
and Morris 39 ultimately conclude that this source Is 
not Important at subsonic speeds, which Is consist- 
ent with the findings of Moore, 33 who studied the 
phenomenon experimentally by artificially exciting a 
jet under conditions that tended to minimize vortex 
pairing. Unlike Klbens, 25 Moore found that the 
broadband noise was usually Increased rather than 
suppressed by the external excitation. He concluded 
that the Instability wave, while not radiating noise 
directly, acted as a conduit through which energy 
could be transferred to the small-scale turbulent 
motion. 

Concluding Remarks 


Q^j = 5 1j Q 0 , (V.J = 1,2.3) (11) 


where Is the Kronecker delta. It follows from 

Eqs. (9) and (10) that 


1 , j=l 




1 - M cos e 

c /c 

0 00 


4 


Q 


2 

0 * 


which more than compensates for the additional two 
Doppler factors In the denominator of Eq. (8) and 
produces a net azimuthal directivity pattern that 
Is given by three Inverse Doppler factors — In 
excellent agreement with the experimentally observed 
one-third octave directivity patterns. 


The Interpretation of this result Is that the 
reduced directivity of the high-frequency sound Is 
due to Interference between the various quadrupole 
components. Note that the mean-square pressure Is 
now the product of an azimuthal directivity factor 
that depends only on e and a circumferential 
directivity factor that depends only on <p. It Is, 


Comparisons of the results of the previous 
section with experiment suggest that the high- 
frequency solution may remain valid at frequencies 
that are low enough to Include the most energetic 
portion of the jet-noise spectrum. This might also 
be true for other more complex turbulent flows. We 
may therefore be able to calculate the sound gener- 
ated In such flows by finding the high-frequency 
solution for a point quadrupole source moving 
through the appropriate mean flow. Durbin 11 
recently obtained such a solution for a completely 
general mean flow. Up to now our remarks have been 
confined to subsonic flows, though much of what has 
been said also applied to supersonic flows. How- 
ever, some Interesting new phenomena also come Into 
play In these flows. 
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Figure 1. - Plate embedded in turbulent airjet. 
Note that the leading edge of the plate is posi- 
tioned in the region of maximum turbulence 
intensity so that its sound field will exceed 
the background jet noise by the maximum 
amount. 
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Figure 2. - Comparison of causal or leading-edge Kutta condition solution with 
data of Olsen (1976) for Uj • 7000 ft s' 1 . 
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